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Abstract 

We propose a device for studying the Fermi-Hubbard model with long-range Coulomb inter- 
actions using an array of quantum dots defined in a semiconductor two-dimensional electron gas 
system. Bands with energies above the lowest energy band are used to form the Hubbard model, 
which allows for an experimentally simpler realization of the device. We find that depending on 
average electron density, the system is well described by a one- or two-band Hubbard model. Our 
device design enables the control of the ratio of the Coulomb interaction to the kinetic energy of 
the electrons independently to the filling of the quantum dots, such that a large portion of the 
Hubbard phase diagram may be probed. Estimates of the Hubbard parameters suggest that a 
metal-Mott insulator quantum phase transition and a d-wave superconducting phase should be 
observable using current fabrication technologies. 
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Along with quantum computers, quantum simulators promise to offer performance ex- 
ceeding what is possible using only classical physics [ll . A quantum simulator is a purpose- 
built device that simulates a particular quantum many-body problem that is intractable on 
a classical computer. Several experiments and theoretical studies have indicated that such 
a device is possible to build [2|, , [qI . Here we show that a similar device can 

be built for the fermionic Hubbard model using a semiconductor quantum dot (QD) array 
system. The two-dimensional Hubbard model is particularly interesting as it is one of the 
central models used to describe strongly-correlated phenomena such as metal-insulator tran- 
sitions magnetism [llj], and high-temperature superconductivity [3]. Despite decades 
of intensive research, a complete solution remains unavailable due to difficulties in numerical 
and analytic methods. 

The type of device we consider is an undoped GaAs/AlGaAs heterostructure, with a 
two-dimensional electron gas system (2DEG) formed at the interface (Fig. [1]). The 2DEG is 
formed by applying a positive voltage to a metallic top gate (the "global gate" (GG)) |l3 |. 
Our choice of an undoped system, as opposed to a modulation doped system, is important 
for a clean realization of the Hubbard model so that impurities in the system are reduced to 
a minimum. In addition to the GG, a 2D mesh gate (MG) is patterned over a large area (e.g. 
30 /im X 30 fim ) of the device. Applying a voltage to the MG induces an in-situ tunable 
periodic lattice potential, such that an array of coupled QDs is created. The MG is separated 
by an electrically insulating layer from the GG, so that the average electron number in the 
QDs can be controlled independently of the inter-QD coupling [13|. The periodicity of the 
MG considered in this paper is A ~ 0.1 /im, which is an experimentally achievable size 
using current lithography techniques. The 2DEG is formed at a relatively shallow position 
relative to the surface (e.g. a depth ci ^ 30 nm). This is advantageous in order to achieve 
sharp QD trapping potentials, as well as reducing the overall Coulomb repulsion between 
electrons in a single QD via electrical screening from the metal gates. Without screening, 

2 fe2 

the Coulomb interaction is typically much larger than the kinetic energy ^ 2m* x'^ 
our semiconductor system (A = 0.1 fim, e = 13eo is the permittivity, m* = 0.067me is the 
effective electron mass in GaAs). The screening effect allows access to an interesting regime 
of the Hubbard model where quantum phase transition (QPT) phenomena are expected to 
occur, with the Coulomb repulsion U and the nearest neighbor hopping t being on the same 
order. 
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One of the key features of our device is that we consider average electron numbers beyond 
the occupation of the first band formed by the periodic potential. Typical high-mobility 
semiconductor samples have electron densities in the region of ~ 10^^ cm~^, corresponding 
to 10 electrons per QD for a periodicity of A = 0.1 j^m. Past studies 9| have assumed 
the occupation of only the lowest energy band of the imposed periodic potential, which 
corresponds to low electron densities where it is difficult to achieve high mobilities due to 
the presence of impurities. Low-density 2DEGs have also been experimentally observed to 
undergo a metal-insulator-like transition as a function of electron density IJ]. In order 
not to mask the effects of the effective Hubbard model due to such low-density effects 
it is advantageous to work in a high electron mobility regime where the system is in an 
unambiguous metallic state. 

Due to the presence of the electrons in the lower energy bands, the Coulomb interaction 
between two electrons in the higher energy bands will experience a screening effect. In the 
standard procedure of obtaining an effective Hubbard model, one usually ignores the presence 



of the core electrons and considers only the outermost filled band of the periodic array 
In order to take into account of this screening effect at a quantum mechanical level, it is 
necessary to incorporate correlations between electrons in the outer orbitals and the core 
electrons. We can take this screening effect into account by the following procedure. First, let 
us model the 2DEG electrons in a periodic potential Vm{x) = Vq [cos(27rs/A) + cos(27r?//A)] 
by a very general multi-Hubbard model 
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where Cjna is the fermion annihilation operator associated with the site j = {jx,jy), band 
n = [rixjUy) (where nx,ny > 1), and spin a. For each band n we may define a Wan- 
nier basis Wn{x), from which we may define the hopping T{n,n',j,j') = j (Pxw^{x — 
Xj)Ho{x)wn'{x - Xj>) and Coulomb U{ji,j2, 33,3^,^1, n2,n-i,n4) = j (Px j (fx'w*^_^{x' - 
Xj^)w^,^{x — Xj.^)Uc{x,x')wn.^{x — Xj.^)wni{x' — Xj^) matrix elements, where the single elec- 
tron Hamiltonian is Hq(x) = — 2^V^ + VA,/(a;). Due to the presence of the metal gates at the 
surface, we use a screened Coulomb interaction Uc{x,x') = e'^fs{x,x')/4:ne\x — x'\, where 



fs{x,x') = l- \x- a;'|/yC 



X — x'y + {y — y'y + {z + z' + 2dy. The form of the screening 



function /^(a;, x') is derived using the method of images by assuming a uniform metal plate 
at the surface. 

Let us now split the Hamiltonian ([T]) into two parts: one corresponding to "on-site" terms 



and the other corresponding to all the remaining terms -ffsitc-sito = H — i/on-site- All terms in 
-f^sitc-sito contain operators connecting two sites j,j' with j' 7^ j. i/on-sito thus describes an 
array of independent QDs, while -ffsitc-sitc contains the interactions between them. In writing 
we have used the identity T{n,n',j,j') = ^^^^ IntbBz <PkEkne''^'^'^'''^i'\ where Ekn 
is the energy dispersion of the nth non-interacting band, and defined tn = T{n^n^ j , j) 
and Uni,n2,n3,n4 = ^(j , j , J , j , rii , 712, rig, 714) . -f^on-sito has precisely the same form as the 
Hamiltonian of an array of isolated QDs [15]. The only difference here is that the single 
particle basis used here is a Wannier basis, instead of the eigenstates of the confining poten- 
tial. For this reason, we henceforth call the many-body system of electrons interacting via 
Hamiltonian ([2]) on a particular site a "Wannier quantum dot" (WQD). In the limit of very 
large barriers between the dots, e„ coincides exactly with the energy levels of an isolated 



In order to analyze the Hamiltonian ([2]), we have performed an exact diagonalization 
study for electron numbers up to = 30. In the inset of Fig. [2] we show the electron 
addition spectrum calculated according to A{N) = Ej^^i + Ej^+i — 2E^ [15], where E]<q 
is the ground state energy of the A^-particle Hamiltonian. We see a peak structure that 
reflects the shell structure of the WQD, with peaks in A{N) occuring at the magic numbers 
of the WQD, i.e. electron numbers corresponding to completely filled shells. The spectrum 
of a single electron WQD is shown in the main plot of Fig. [2] to show the shell structure. 
Unlike normal QDs, the shells are typically only doubly degenerate at most. The two-fold 
degeneracy originates from the x-y symmetry of the potential. For electron numbers in the 
WQD corresponding to the doubly degenerate shells, we have verified that Hund's rule holds 
according to our numerical calculations, with ground states occuring for a total 2;-component 
of spin Sz = ±1, rather than 5"^ = at half-filled shell fillings. 

Using the diagonalized states of ([2]), we may write down an effective Hubbard Hamilto- 
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nian. The approximation we make is that for a given electron density, only one shell of a 
WQD is relevant to describe the low energy physics. For fillings corresponding to a non- 
degenerate shell, we may make the state associations \0)j = \Nb, 0, j), | t)j = \Nb + l, 1/2, j), 
I l)j = \Nh + 1,— l/2,j), I t|)j = lA^";, + 2,0,j) for each site j, where the eigenstates of 
on a single site j are denoted |A^, S^,j), and the total number of electrons in the lower 
energy "core" shells is called the base electron number A*";,. For doubly degenerate shells, 
a maximum of four electrons may occupy the shell, which gives rise to a two-band Hub- 
bard model. In this paper, for simplicity we consider the non-degenerate case, although the 
generalization to the two-band case may be performed straightforwardly. 

Let us now introduce a particle number operator nj^ acting on the states Sj = {\0)j, \ f 
I I H)j} with the properties nj^\0)jf = 0, nj^\a')j> = djfdcra'W), Ti)i' = Sjj'l Ti 
). In the space of Sj, we may write down an effective Hamiltonian 

H:L.. = J2^Eo+ Y1 f^aUj^ + Unj^nj^) (3) 

j <T=±l/2 

where Eq = Ejy^fi, fj^o- = -EAr^+i,cr — -Eat^^o, U = Eiy^+2,o — J2a=±i/2 En^+i,^+ ^N^fi- The /io- is an 
effective chemical potential term, while U is an effective Coulomb on-site repulsion energy. 
Figure [3] shows the on-site repulsion as a function of the periodic potential amplitude for 
a base electron filling A*";, = 0. We see that the periodic potential amplitude increases the 
effective on-site interaction, in agreement with previous calculations jol. 

The hopping terms between sites are calculated according to the transition that the inter- 
site Hamiltonian i^sitc-sito induces between the WQD eigenstates by making a transformation 
of the hopping terms in H^i^^,^^^^ into the truncated basis states Sj. To a good approximation 
we find that the hopping can be approximated by the hopping integral of the band corre- 
sponding to the outermost shell of the non- interacting WQD. The nearest neighbor hopping 
amplitudes for various base fillings Nh are plotted in Fig. [31 We see that as the periodic po- 
tential is increased, the hopping is suppressed, due to the increased strength of the barriers 
between the QDs. As the number of electrons is increased, the tunneling is enhanced due 
to the electrons occupying higher energy bands. 

Table [T] gives our calculated values of various Hubbard parameters for two potential 
amplitudes Vq = 0.56 and 5.4 meV and various base electron numbers Ni,. For a given 
potential Vq, the nearest neighbor and next-nearest neighbor hopping amplitudes t and t' 
increase with the base filling A";,. Naively one may conclude based on this that using high 
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base electron numbers is always advantageous due to the larger energy scale of the effective 
Hamiltonian. However, one must also consider whether inter-band transitions are suppressed 
sufficiently in order that the single band (or two-band) approximation is valid. For low 
temperatures k^T <^ 2m* x'^ ' approximate criterion is when there is no energy overlap 
between a given energy band and all other bands for the non-interacting band spectrum. 
For Vq = 0.56 meV, a simple calculation reveals that only the lowest band is separated in 
energy from the other bands, hence an effective Hubbard model will only be formed for 
A*";, = 0. For Vq = 5.4 meV, we find that the four lowest energy bands are separated, hence 
we require < 12. Since both Vq and the average electron density are freely choosable 
parameters, in practice we may always obtain an effective Hubbard model by appropriately 
choosing these values. 

By choosing a band with a small U/t and increasing the periodic potential Vq, we expect a 
metal- Mott insulator transition to occur in the system ^, [lo|. In our proposed device, such 
a transition may be identified by measuring the zero-bias differential conductance across the 
source and drain contacts. The conductance as a function of the chemical potential has been 
theoretically investigated in previous studies [16|]. In the metallic state, we expect the system 
to be conducting when the Fermi level lies within a band created by the effective Hubbard 
model. In the Mott insulating limit, the spectrum should reduce to that of a large array of 
isolated QDs. The characteristic conductance spectra of the metallic and insulating states 
differ in each case, thus by varying the potential Vq at a fixed Fermi level, one should be 
able to distinguish a transition between the two phases. Magnetoconductance measurements 
provide an independent check of a QPT between the metallic and Mott insulating states. 
The phase-breaking length will reduce to a length comparable to the QD size in the case 
of a Mott-insulating QD array [fT^. Experimental data consistent with these expectations 
have already been observed to some extent 
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19j. However, a clear metal-Mott insulator 



transition has never been identified to our knowledge. 

In order that such QPT phenomena are not washed out due to the temperature effects, we 
require the Hubbard parameters to be larger than the thermal energy f/,t > fcijT l2Qi. Using 
dilution refrigerator techniques, temperatures of T = 10 mK are reachable, corresponding 
to a thermal energy of ksT ^ 1 /leV. Previous theoretical estimates have suggested that 
temperatures of T ^ O.lt are necessary to observe an antiferromagnetic (AF) phase, and T ^ 
0.02t for a rf-wave 
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we require a hopping of about t > 0.01 meV for observation of an AF phase and t > 0.04 meV 
for observation of a SC phase. Comparison of these numbers with Tabled we find that both 
these phases should be possible to observe within the Hubbard approximated bands. The 
AF nature of the insulating phase can be determined from temperature-varying magnetic 
susceptibility measurements 22] . Evidence of Cooper-pair formation may be obtained from 
the magnetocapacitance oscillation period, by observation of the Cooper pair charge of 2e 
in the strongly coupled QD array. 

In summary, we have proposed an experimentally viable quantum simulator for the one- 
and two-band Hubbard models using a semiconductor QD array device. For a given average 
electron number in the QDs, the low-energy physics may be described by an effective one- or 
two-band Hubbard model. Our scheme may be easily generalized to different lattice geome- 
tries simply by adjusting the mesh gate design and voltage. Examining the region U/t ^ 1 
produces an effective t-J or Heisenberg model, while spin models involving frustration may 
also be explored by fabricating a triangular lattice mesh gate. Another possibility is to 
introduce controlled disorder into the system by randomly varying the mesh dimensions in 
the lattice, and thereby producing a Hubbard-Anderson model. 

This work is supported by JST/SORST, NTT, the University of Tokyo, and the Special 
Coordination Funds for Promoting Science and Technology. We would like to thank M. 
Beasley, I. Fisher, M. Jura, M. Topinka, J. Berengut, and P. Recher for helpful discussions. 
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TABLE I: (a) The on-site Coulomb energy U, the nearest-neighbor Coulomb energy V, nearest- 
neighbor hopping t, and next nearest-neighbor hopping t' are computed for various base filling 
numbers A^^ for a potential amplitude Vq = 0.56 meV. The nearest neighbor Coulomb interaction 
V is calculated using the standard single band approximation. Numbers in brackets denote errors 
on the last significant digit. The band index n corresponding to the filling level of the non- 
interacting Wannier quantum dot is given. For base fillings corresponding to a two-band Hubbard 
model, we take the average on-site Coulomb interaction between the electrons species, (b) As for 
(a) but for Vq = 5.4 meV. 

Nb n U (meV) V (meV) t (meV) t' (meV) 



(1,1) 


0.95(3) 


0.24(7) 


0.080 


0.011 


2 (1,2) (2,1) 


0.8(3) 


0.21(4) 


0.192 


0.028 


6 (2,2) 


1.5(10) 


0.14(3) 


0.305 


0.044 


8 (1,3) (3,1) 


1.3(10) 


0.24(3) 


0.323 


0.021 


12 (2,3) (3,2) 


1.3(5) 


0.5(4) 


0.435 


0.037 


16 (1,4) (4,1) 


1.0(5) 


0.6(4) 


0.438 


0.020 



Nb n U (meV) V (meV) t (meV) t' (meV) 



(1,1) 


4.3(1) 


0.17(7) 


0.0016 


2.2 X 10-6 


2 (1,2) (2,1) 


3.2(2) 


0.10(3) 


0.020 


4.8 X 10-^ 


6 (1,3) (3,1) 


2(1) 


0.30(6) 


0.129 


0.016 


10 (2,2) 


2(1) 


0.05(3) 


0.038 


9.7 X 10"^ 


12 (1,4) (4,1) 


4(2) 


0.3(1) 


0.32 


0.035 


16 (2,3) (3,2) 


2(2) 


0.10(3) 


0.15 


0.017 
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FIG. 1: (color online) (a) A schematic top view of the proposed device. A two-dimensional Schottky 
mesh gate (MG) is patterned in the central region, and a top global gate (GG) is placed on top 
of the MG separated by an insulating layer such as Si3N4 (INS). Source (S) and drain (D) ohmic 
contacts in the Hall bar mesa structure access the two-dimensional electron gas (2DEG) system, 
(b) Cross-sectional view of the device. The 2DEG is formed at the interface between AlGaAs and 
GaAs by applying a voltage to the GG. The periodic potential is created by the 2D MG. 
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FIG. 2: The single particle energy spectrum of a 2D Wannier quantum dot (WQD) versus the 
energy level i. The ith level can be occupied by two electrons due to spin. The inset shows the 
addition energy A{N) of a WQD. We use parameters A = 0.1 /xm, Vq = 1.1 meV, and d = 10 nm 
for the calculation. 
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FIG. 3: The on-site Coulomb energy U (dashed line, right axis) for Nj, = and hopping t (solid 
lines, left axis) versus the periodic potential amplitude. The base electron number for each t are 
labeled. Parameters A = 0.1 //m and d = 30 nm are used for the calculation. 
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